Formation of singularities on the surface of a liquid metal in a 

strong electric field 

N. M. Zubare 

Institute of Electrophysics, Ural Branch, Russian Academy of Sciences, 
106 Amundsen Street, 620016 Ekaterinburg, Russia 

The nonlinear dynamics of the free surface of an ideal conducting liquid in a strong 
external electric field is studied. It is establish that the equations of motion for such a 
liquid can be solved in the approximation in which the surface deviates from a plane 
by small angles. This makes it possible to show that on an initially smooth surface 
for almost any initial conditions points with an infinite curvature corresponding to 
branch points of the root type can form in a finite time. 



I. INTRODUCTION 

A flat surface of a conducting liquid placed in a strong external electric field is known 
to become unstable if the field strength E exceeds a critical value E c 2 = 8iTy/gap, 
where g is the acceleration of free fall, a is the surface tension, and p is the density of the 
medium. The interaction of the electric field and the charges induced by this field on the 
surface of the liquid causes surface perturbations to grow rapidly and regions of substantial 
curvature to form in a finite time jj, 0|. The result may be an increase in the field energy 
density at the surface, initialization of emission processes, and, finally, vacuum breakdown 
fl. Moreover, there are indication that the liquid phase p.a y8 an important ro.e in the initial 
stages of explosive electron emission p|. All this means that one must build a meaningful 
theoretical model of the nonlinear stages in the development of an instability, a model that 
can be used to describe the formation of a singular profile of the surface of the medium (a 
liquid metal in applications). 

The present paper studies the nonlinear dynamics of an electrohydrodynamic instability 
in the limit of a strong electric field, E 3> E c , when both surface tension and gravity 
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can be ignored. The interest in this limit is due, in particular, to the recent discovery 
of systems with anomalously low critical fields, E c ~ lkVcm -1 (Ref. ^). The nonlinear 
stages in the development of an instability are studied by perturbation techniques that use 
series expansions in a small parameter, the angle of the slope of the surface. Of course, 
the introduction of such a small parameter makes it impossible to describe the formation 
(observed in experiments) of strong singularities, with corresponding slope angles of order 
unity. Nevertheless, using the model adopted in this paper, we can show that for almost 
any initial conditions at the surface of the conducting liquid, it takes only a finite time for 
points with infinite curvature to form on the surface. Thus, even in the weakly nonlinear 
stages in the development of a nonlinear ity there is the tendency for a singular profile of the 
liquid surface to form. 

The plan of the paper is as follows. In Sec. II we derive the main equations of motion, 
which describe the evolution of the free surface of an ideal conducting liquid in a strong 
external electric field. In Sec. Ill we use the approximation of small angles characterizing the 
slope of the surface to build a nonlinear model of the development of an electrohydrodynamic 
instability. Section IV is devoted to a study of the dynamics of one-dimensional surface 
perturbations. Integration of the model equations shows that it takes only a finite time for 
weak singularities of the root type to form in the system, i.e., singular points at which the 
curvature of surface is infinite (see also the Letter [9|). 

Mathematically, the formation of singularities can be explained by the violation of the 
analyticity of the complex velocity potential due to the movement of singularities, or branch 
points, to the boundary. On the whole, such behavior is similar to that of the velocity 
potential of an ideal liquid in the absence of external forces jlO, 11, In Sec. V we use the 
example of the evolution of single perturbations to show that the formation of singularities 
occurs before the small-angle condition is violated because of the development of a linear 
instability (the branch point of the root type agrees with the small-angle approximation). 
In Sec. VI we study the behavior of the boundary of a liquid metal under the assumption 
that self-similarity is retained in a small neighborhood of a singularity in the crossover from 
one-dimensional perturbations of the surface to arbitrary perturbations. Finally, in Sec. VII 
we discuss the role that branch points of the root type play in the evolution of the system. 
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II. INITIAL EQUATIONS 

Consider the motion of a conducting liquid that occupies the region — oo < z < r)(x, y, t) 
and is subjected to a strong electric field E. We assume that this liquid is ideal and its 
motion is vortex-free. Then the potential $ of the liquid velocity is determined by the 
time-dependent Bernoulli equation 

where p is the pressure and F is a function of time. Moreover, for potential flow of an 
incompressible fluid we have A$ = 0. The equations of motion must be augmented by the 
kinematic condition at the free surface, 

th= [$*-V77-V$]| z=r?) 

by the condition at infinity, V$|^_ 00 —>■ 0, and the condition imposed on the pressure at 
the conductor- vacuum boundary [3], 

[87rp+(Vy?) 2 ]| 2=r? = 0, 

where ip is the potential of the electric field. 

The electric potential in the absence of space charges is described by the Laplace equation 
Aip = together with the conditions that everywhere on the surface of the conductor the 
potential be the same, ip\ z= „ = 0, and that the field be uniform at infinity, ^L-^oo — > —Ez. 

Note that these equations of motion have a Hamiltonian structure and the functions 

,(«,»,*) and « x ,y,t) = »U - canonicaUy conjugate HQ: 

dip 8H dr] _ 5H 
dt Si] ' dt Sip' 

where the Hamiltonian 



2 J 8vrp 

z<r) z~>7] 

coincides, to within a constant, with the total energy of the system. 



III. THE SMALL-ANGLE APPROXIMATION 



Using Green's formulas, we can write the Hamiltonian in the form of the surface integral: 



H 



> d<5> E 2 r] dip' 
2 dn 87rp dn 



ds, 
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where <p = (ip + Ez)/E is the perturbation of the scalar potential, ds is the surface area 
element, and d/dn is the normal derivative at the surface s. 

From now on we assume | <C 1, which corresponds to small surface-slope angles. This 
allows expanding the normal derivatives in powers series of the canonical variables. Then 
for the Hamiltonian we have 



H=J^ (f+kf- 1 ^ - V?7 • f+VT '- V) d 2 r 
-J {fAfl 1 ^ + Vr] • f^Vfr 1 ^ d 2 r. 



Here k is the two-dimensional integral operator with a kernel whose Fourier transform is 
equal to the absolute value of the wave vector: 

+00 +00 



kf = - 



1 



f(x',y') dx'dy' 



2W J [( x /_ x )2 + ( y /_ y )2]3/2' 

00 —00 



The nonlinear operators T± defined as 



{±r]) n k r 



n=0 



act as shift operators (i.e., f\ z=v = Tf\ z= o) for harmonic functions that decay as z — > =Foo. 
If we limit ourselves to second- or third-order terms and introduce scaling 



t -> tE- L {A7Tp) 1/2 , i{j -> ^/(4vrp) 1/2 , # -> HE 2 /(4tt P ), 



1/2 



we arrive at an expression for the Hamiltonian: 

H=^J ipkip -r]kr] + r] ((W>) 2 - (^) 2 + (V??) 2 - (fc^) 2 ) 
The equations of motion corresponding to this Hamiltonian are 



d 2 r. 



g t + kg = l [{kff ~ (V/) 2 + (kg) 2 - (Vg) 2 ] + k (/ - g)kf + V • [(/ - <?)V/] , 



(1) 



(2) 



+ V • [(g - f) Vg] , (3) 



where we have changed from the variables 77 and ^ to the normal variables / and g: 



f = 



1p + Tj 



9 = 



lp — Tj 
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In the linear approximation, Eq. (2) describes the relaxation of g to zero with a char- 
acteristic times l/\k\. In this case in the right-hand sides of Eqs. (2) and (3) we can put 
g — 0, which means we are examining the perturbation-buildup branch with allowance for 
a quadratic nonlinearity. This leads us to the following system of equations: 

9t + kg = \ (kff ~ \ (V/) 2 + k(fkf) + V • (/V/), (4) 

ft ~ kf = \ (kf f - \ (V/) 2 . (5) 

Thus, we have shown that studying the dynamics of perturbations of the surface of a 
conducting medium in a strong electric field in the small-angle approximation amounts to 
studying the system of equations (4) and (5). What is important about this system is that 
the nonlinear equation (5) does not contain the function g and that Eq. (4) is linear in g 
and can easily be solved: 

t +oo +oo 

= J_fff G(x',y',t') (t-t')dx'dy'dt' 

9 2ttJ J J [( x /_ ;r )2 + ( y /_ y )2 + ( f _^2]3/2' 

— oo — oo 

G(x, V,t)= l - {kff - l - (V/) 2 + k(fkf) + V • (/V/), (7) 
where we assumed that <?|t=o = 0. 



IV. FORMATION OF A BRANCH POINT OF THE ROOT TYPE 

In the case of one-dimensional perturbations of the surface (we ignore the dependence 
of all the quantities on y), the integral operator k can be expressed in terms of the Hilbert 
operator H: 

+oo 

i = -» A, ^I/i<^. 

OX IT J X' — X 

— oo 

Then the model equations (4) and (5) can written 

9t ~Hg x = \ (Hf x f - l - (f x f + H(fHf x ) x + (ff x ) x , (8) 



/* + Hf* = l (Hf x f - \ (f x f. 



(9) 
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For further discussions it is convenient to introduce functions that are analytic in the 
upper half-plane of the complex variable x: 

= Pf, x = Pg, 

where P — (1 — iH)/2. Since applying the Hilbert operator to a function that is analytic in 
the upper half-plane amounts to multiplying that function by the unit imaginary number, 
Eqs. (8) and (9) take the form 



Re (<f> t + i<p x + <Pl) = 0, 
Re (xt - iX, + &-2P (00*0,) = °- 



Thus, the integro-differential equations (8) and (9) can be studied simply by analyzing the 
inhomogeneous linear equation 

X*-% = -0'+2P(00 :E ) x (10) 
and (separately) the nonlinear partial differential equation 

</h + i<f>x = -<f>l- (11) 

For the sake of convenience we introduce a new function, v = <p x . In terms of this function, 
Eq. (11) becomes 

Vt + iv x = —2vv x . 



Note that this equation coincides with the one proposed by Zhdanov and Trubnikov [15|, 1||, 
who used it to describe the nonlinear stages in the development of tangential discontinuities 
in hydrodynamics. More than that, if we replace x by x — > x + it, we arrive at the equation 



derived in Refs. 110, IllL Il2« as a result of a discussion of the nonlinear dynamics of a free 
surface of ideal liquid in the absence of external forces, where it describes the temporal 
evolution of the complex-valued velocity. The solution of this first-order partial differential 
equation can be found by using the method of characteristics: 

v = Q(x'), (12) 



x = x' + it + 2Q(x')t, 



(13) 
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where the function Q is determined by the initial conditions Q(x) = v\ t=0 . 



Let us show, by analogy with Refs. [1C , 
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121 ] . that these relations describe (if we require 



that v be analytic) the formation of a singularity in a finite time. The problem of finding 
the explicit solution reduces to analyzing the map x — > x' specified by Eq. (13). Generally, 
this map ceases to be one-to-one at points where 



dx/dx' = 1 + 2Q x ,t = 0. 



(14) 



The relationship (14) specifies a path x' = x'(t) in the complex x' plane. Then the motion 
of the branch point of the function v is given by 

x(t) = x'(t) +it + 2Q(x'(t))t. 

At the time to when the branch point reaches the real axis the analyticity of v is violated 
and the solutions of Eq. (9) become singular. 

Let us examine the behavior of the solutions near a singularity. Expanding (12) and (13) 
in a small neighborhood of the point t = t , x = Xq = x(to), (t ), in the leading 

order we get 

v = qo — Sx' /(2to), 5x = iSt + 2q 5t + q"to(Sx') 2 , 

where q = Q(x' ), q" = Q x ' x '(x' ), 5x = x — x , 5x f — x' — x' , and 5t = t — to- 
Excluding Sx' from these expressions, we obtain 

5x — i5t — 2q 5t 1 ' 



Aq"tl 



v = q - 

The derivative of this expressions with respect to x is 

A ~ - - - [I6q%(8x - iSt - 2q 5ty ~ 1/2 



(15) 



(16) 



which shows that cj) xx behaves in a self-similar manner (Sx ~ 5t) and becomes infinite as 
5t -> 0. 

As for the complex- valued function x, the equation that determines its temporal dynamics 
(Eq. (10)) can be integrated by the method of characteristics (see Eqs. (6) and (7)). Taking 
the initial condition in the form x\t=o = yields 
t 



\ J D(x + it — it' ,t') dt' ', D(x 1 t) 
o 



+ 2P 
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Inserting (15) into this expression, we see than near the singularity the derivative Xxx can 
be expressed in terms of <p xx : 

= (go ~ go) , 

Xxx i •% Yxx- 

(go + l) 

This means that the analyticity of Xxx is violated at time t = t . 

How does the surface of the liquid metal behave at the time when the singularities develop 
in the solutions of Eqs. (10) and (11)? Allowing for the fact that rj = f — g, we find that 
the surface curvature 

K = r) xx (l + r)l)- 3/2 
is specified, to within a quadratic nonlinearity, by the expression 



K ^r] xx = 2Re (4> xx - x 



XX I 



Substituting the expression for <f> xx and Xxx found earlier, we find that in a small neighbor- 
hood of the singular point 

(go - go 



K « 2Re 
Since 4> xx is given by (16), we have 



(go + i) 



(17) 



K\ X=X0 ~ \5t\- 1 ' 2 , K\ t=t0 ~ \5x\' l/ \ 

i.e., it takes a finite time for a singularity of root type (branch point) to form at the surface, 
and the curvature of the surface of the liquid at this point is infinite. 

To conclude this section we note that since we have ip — f + g, near the singularity a 
relationship holds for the complex- valued potential of the liquid flow, \P = 2Pip: 

, (go + i) , 



^xx = 2{4>xx + Xxx) ~ 2- 



■'xx 1 



(go + i) 

i.e., the first derivative of the complex- valued velocity also exhibits singular behavior as 
5t — > 0. This means that, as in Refs. 10, 11, the formation of singularities can be 
interpreted as the result of violation of the analyticity of the complex-valued potential due 
to the movement of the singularities of the potential to the boundary. 



V. EVOLUTION OF A SINGLE PERTURBATION 



We use a simple example to show that at the time when a singularity in the solutions of 
Eqs. (8) and (9) develops the applicability conditions for our model are met. 



9 



We take the initial perturbation in the form 

f\ t=0 = -ea m k m - l \n(x 2 + a 2 ), (18) 

where m is a positive integer, and the parameters a and e take positive values a > and 
e > 0. This situation corresponds to a one-dimensional single perturbation of the surface 
symmetric with respect to point x — 0, at which the surface curvature is negative. The 
characteristic slope angles of the surface are determined by the parameter e, which we 
assume small. 

Note that in the linear approximation Eq. (9) becomes 

ft + Hf x = 0. 

Its solution with the initial condition (18) is 

f(x, t) = -ea m k m - 1 In {x 2 + (a - t) 2 ), 

i.e., within the linear model the perturbation grows without limit and becomes infinite at 
the time t = a, which of course violates the applicability conditions for this model. 

Will introducing nonlinearity into the model permit a singularity to develop in the so- 
lution before the condition \r/ x \ ps \f x \ <C 1 breaks down? (The branch-point nature of this 
singularity agrees with the condition that the angles be small.) To answer this question, we 
will examine the evolution of the perturbation (18) according to the nonlinear equation (9). 

The symmetry of (18) implies that the singularity forms at point x — 0. Then from (13) 
and (14) it follows that the time t at which the singularity develops can be found by solving 
the following equations simultaneously: 

x' + it + 2Q(x' ) t = 0, 1 + 2Q x ,{x' ) t = 0, 

where the function Q corresponding to (18) has a pole of order m at the point x' = —ia: 

x 1 + la J 

Expanding in a power series in the small parameter e, we obtain to leading order the fol- 
lowing: 



to ~ a 



m+1 ,xl/(m+l) 

1 (2eml) M ' 

m 



Xrt ~ — Id 



1 - (2em\) 



l/(m+l) 
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Since in the linear approximation the singularity is formed at time t = a, the above expres- 
sion for t implies that the nonlinearity accelerates the formation of the singularity (but if 
e < holds, the nonlinearity delays the onset of the instability). 

Plugging the above expression for x' Q into the expression for Q and its second derivative 
Qx'x', we obtain 

q» „ -^±±1 (2em!)- 1/(m+1) , q « (25m!) 1/(m+1) . 

Thus, for perturbations of the form (18) the parameter q" is finite. This means that the 
dynamics of surface perturbations near a singular point is described fairly accurately by Eqs. 
(15)— (17). As for the parameter q , the smallness of e implies \q \ <C I. This is an important 
result. The important point is that this parameter, as (12) and (15) imply, determines the 
characteristic angles of slope of the surface by the moment of singularity formation. Then 
for the derivative r\ x at the time of collapse the following estimate holds: 

~ e 1 ^ « 1, 

i.e., the characteristic angles remain small, even through they increased by a factor of 
£ -m/(m+i)_ j n ^Yas case there is not enough time for the small-angle condition to be vio- 
lated as a result of the development of a linear instability, and the proposed model (Eqs. 
(8) and (9)) is closed in the sense that if the initial perturbation meets all the conditions 
needed for the model to be valid, this property is retained throughout the entire evolution 
until the time of collapse, t . 

We now discuss the behavior of a perturbation of the electric field at the conducting 
surface, 

5E(x,t) = -E-^- 

9n z=v 

near the singularity. Clearly, in the linear approximation the field perturbation is specified 
by the formula 

5E « —E Hf x = 2Elm(v). 
Substituting v of Eq. (15) in this expression, we find that at the singular point 

5E\ Sx =8t=o ~ 2Elm(q ). 

Since the parameter q is small, the perturbation of the electric field at the time of singularity 
formation remain much weaker than the external field (both 5E X and SE t are singular). 
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VI. SELF-SIMILAR SOLUTIONS IN THE GENERAL CASE 



In Sec. 5 we found that \q \ <C 1 holds in the small- angle approximation. This allows 
us to ignore the dependence on this parameter in Eqs. (16) and (17). Then, if q" ^ 0, the 
dynamics of the surface curvature in the one-dimensional case is given by the expression 

K « -Re [Aq"tl(5x - iSt)] ~ 1/2 . (19) 

In particular, for a perturbation that is symmetric with about the point x = xq we have 

1/2 



K 



-St + V5x 2 + 5t 2 
5x 2 + 5t 2 



Here the initial conditions determine only the constant factor. Thus, the behavior of the 
system near a singular point is of a universal nature. 

Let us take a particular case q" = 0. Suppose, for instance, that 



d 2 Q 



dx' 2 



c'—x' 



dx 



tn—l 



o, 



9n 



c'—Xq 



dx 



in 



c'—x' 



where n > 2. Then, expanding (13) in a power series in 5x' up to the nth power, in the 
leading order we get 



1 ^ / n\ 
K ~ Re 

nt \2t q n 



l/n 



(8x - iSt) 1 '"' 1 . 



(20) 



The formulas (19) and (20) show that for an arbitrary one-dimensional perturbation of 
the surface satisfying the condition j^l <^ 1 the curvature near the singularity behaves 
self-similarly: 

K^J-hf^j), (21) 
\5t\p \\6t\J' v ; 

where h is an unknown function, and the exponent is given by the expression 



p = (n — l)/n. 



(22) 



with n a positive integer. 

Note that (19) and (20) are the exact solutions of the linear integro-differential equation 

K t + HK X = 0, 

which describes the evolution of the surface curvature in the linear approximation, as follows 
from (9) with allowance for the fact that K = f xx holds in the leading order. For an arbitrary 
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exponent p, Eq. (21) specifies the class of self-similar solutions of this equation. This means 
that, on the one hand, the dynamics of the surface near a singularity is described by the 
self-similar solutions of the linearized equations of the model and, on the other hand, that 
the presence of a nonlinearity leads to a situation in which of all the possible self-similar 
solutions only those with rational values of p specified by the condition (22) are realized 
(from general considerations it follows that p = 1/2). 

It is therefore natural to assume that in the two-dimensional case, as in the one- 
dimensional, the solutions in a small neighborhood of the singularity are self-similar: 



where p satisfies the condition (22). A characteristic feature of the weak-nonlinearity ap- 
proximation in our problem is that the specific form of the dependence of all quantities on 
the self-similar variables can be treated using the equation 



whose linearity makes it possible to effectively study the formation of two-dimensional sin- 
gularities at the surface of a conducting medium. Note that this statement is valid if p < 1, 
which, as condition (22) shows, is met in our case in a natural manner. The point is that at 
p = 1 an expression of the form (23) corresponds to the symmetries of the initial nonlinear 
equation of motion. This mean that near a singularity the contribution of a nonlinearity is 
comparable to that of the linear terms, and the analysis of the behavior of the surface lies 
outside the scope of this paper. 

Substituting (23) in (24), we arrive at the following integro-differential equation for the 
unknown function h: 



where £ = Sx/\5t\ and ( = 5y/\5t\. Since the profile of the surface begin to form at the 
periphery and only then is propagated to the point 5x = 5y = 0, at the time of collapse 
the curvature of the surface in a small neighborhood of the singular point is determined by 
the asymptotic solutions of this equation as £ 2 + ( 2 — > oo. As can easily be shown, these 
solutions are described by the partial differential equation 




(23) 



K t = kK, 



(24) 



Sh + {h c + P h = k(s,{)h, 



^ + (h c + ph = 0, 
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whose general equation is 

h=[e+c 2 ]- p/2 nc/o, 

where F is an unknown function. Plugging this expression into (23) and introducing polar 
coordinates, 

Sx = r cos /3, 5y = r sin (3, 

we arrive at the following formula for the curvature of the surface near the singular point: 

F(tg/3) 



K 



t=t ~ 



We see that we are again dealing with a branch point of the root type. 

VII. CONCLUSION 

Our analysis of the evolution of the boundary of a conducting liquid in a strong electric 
field within the small-angle approximation has shown that for almost any initial conditions 
on an initially smooth surface the presence of a nonlinearity gives rise to points at which 
the curvature of the surface becomes infinite. These points correspond to branch points 
of the root type. However, the presence of such singularities does not ensure a significant 
concentration of the electric field near the surface of the conductor and, hence, cannot 
by itself lead to vacuum breakdown. In this case we may assume that the main role of 
these branch points in the general evolution of the system is, in time, to generate stronger 
singularities capable of substantially influencing the emission from liquid metal; in particular, 
capable of ensuring the conditions needing for the initiation of explosive electron emission. 
Among such singularities are, for instance, discontinuities in the first derivative of the surface 
profile, which were observed in experiments Note that a theoretical study of such 

singularities lies outside the scope of the small-angle approximation and requires allowing 
for surface tension. Indeed, the applicability of the model adopted in the present paper 
is limited to the scales \5x\ ^> I, where the parameter I is the characteristics length on 
which the capillary effects become important. This parameter, on dimensional grounds, is 
determined by the ratio of surface tension to the electric- field energy density: I ~ a/E 2 . 
The stabilizing effect of the surface pressure means that at time to the curvature of the 
surface is still finite (K ~ 1 //) and the profile of the surface is smooth, and the formation 
of a singular profile begins in the later stages in the development of an instability. 
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